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ABSTRACT 

We study the star formation rates (SFRs) of galaxies as a function of local galaxy density at 0.6 < 
z < 0.9. We used a low-dispersion prism in IMACS on the 6.5 m Baade (Magellan I) telescope to obtain 
spectra and measured redshifts to a precision of (7^/(1 + z) ~ 1% for galaxies with zab < 23.3 mag. We 
utilized a stellar mass-limited sample of 977 galaxies above M > 1.8 x 10 10 M© (log M/Mq > 10.25) 
to conduct our main analysis. With three different SFR indicators, (1) Spitzer MIPS 24 /im imaging, 
(2) SED fitting, and (3) [OII]A3727 A emission, we find the median specific SFR (SSFR) and SFR to 
decline from the low-density field to the cores of groups and a rich cluster. For the SED and [Oil] 
based SFRs, the decline in SSFR is roughly an order of magnitude while for the MIPS based SFRs, the 
decline is a factor of ~ 4. We find approximately the same magnitude of decline in SSFR even after 
removing the sample of galaxies near the cluster. Galaxies in groups and a cluster at these redshifts 
therefore have lower star formation (SF) activity than galaxies in the field, as is the case at z ~ 0. 
We investigated whether the decline in SFR with increasing density is caused by a change in the 
proportion of quiescent and star forming galaxies (SFGs) or by a decline in the SFRs of SFGs. Using 
the rest-frame U — V and V — J colors to distinguish quiescent galaxies from SFGs (including both 
unattenuated blue galaxies and reddened ones) we find the fraction of quiescent galaxies increases 
from ~ 32% to 79% from low to high density. In addition, we find the SSFRs of SFGs, selected based 
on U — V and V — J colors, to decline with increasing density by factors of ~ 5 — 6 for the SED and 
[Oil] based SFRs. The MIPS based SSFRs for SFGs decline with a shallower slope. The declining 
SFRs of SFGs with density is paralleled by a decline in the median A v , providing indirect evidence 
that the cold gas content that fuels future SF is diminished in higher density environments. The order 
of magnitude decline in the SSFR-density relation at 0.6 < z < 0.9 is therefore driven by both a 
combination of declining SFRs of SFGs as well as a changing mix of SFGs and quiescent galaxies. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: clusters: general — galaxies: 
clusters: individual (RX J0152.7-1357) 



1. INTRODUCTION 

Hierarchical structure formation results in more mas- 
sive dark matter halos assembling at late times. Simu- 
lations predict that halos that host L* galaxies, groups, 
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and clusters, grow in mass by fa ctors of ~ 2 — 3 be- 
tween z ~ 1 and z ~ (e.g., iWechsler et al.l 120021 : 
Ivan den TXo sch 2002). At z ~ 0, these massive halos, or 
high density environments, typically host more evolved 
galaxy populations. The morphological study of lDresslerl 
(1980) provided the first quantitative glimpse of this 
fact, finding elliptical and SO, or early type galaxies, to 
dominate in high-density regions of clusters while spirals 
and irregulars, or late-types, dominate at lower densi- 
ties. This morphology density relation (MDR) has also 
been found to extend to group environments at z ~ 
([Postman fc Gcllcr Il984h . These high density environ- 
ments have continuously accreted galaxies from a range 
of lower density environments over time. Thus, mass 
growth in such high density environments would lead 
galaxies within newly accreted halos at z ~ 1 to undergo 
various transformations in their properties if the general 
environmental trends at z ~ hold out to these higher 
redshifts. 

In addition to morphology other galaxy properties 
such as rest-frame colors (e.g., iKauffmann et all 1200 
iHogg et al.|[200llBlanton et al.ll2005HBaldrv et al.ll200 
and st ar formation rate s (SFR s, e.g., Abrahar net all 
19961: iHashimoto et all 119981 : iGomez et all 120031 : 
Kauffman n et all [2004) have also been found to cor- 
relate with the local galaxy density at z ~ 0. Large 
rcdshift surveys find such correlations in the field as 
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well, indicating the importance of mechanisms that 
operate at halo mass scales below that of clusters in 
transforming galaxy properties at z ~ 0. In general 
the works above find red, quiescent, early-type galaxies 
represent a larger proportion of the population in higher 
density regions at z ~ such as the cores of groups 
and clusters, while blue, star-forming, late-type galaxies 
make up a larger share in lower density regions. 

While several studies at z ~ have found the mean 
or median SFRs of galaxies to be lowe r at higher densi- 
ties s u ch as in groups and cl usters (e.g., [H ashimot o • et al.l 
[1991 iGomez et all l200l iKauffmann et all l2004D ~at 
higher redshifts the direction of this SFR-density trend 
has been the subject o f deb ate. For luminosity- l imited 
samples, lElbaz et al.l (|2007t ) and iCooper et atl (|2008l ) 
find that at z ~ 1 the SFR-density relation reverses, 
such that galaxies at higher densities have higher SFRs. 
Meanwhile, for a stellar mass-limited sample at z ~ 0.8, 
iPatel et al.l (|2009bf ) found the obscured SFR-density re- 
lation to decline at higher densities, much like at z ~ 0. 
The direction of the SFR-density trend at z ~ 1 has 
implications for the effectiveness of the various physical 
processes that operate in different environments in reg- 
ulating star formation (SF) and evolving galaxies. For 
example, a reversed SFR-density trend at z ~ 1 would 
imply that the physical processes that lead to the shut- 
down of SF in high density group environments produce 
fewer red, non-star forming galaxies relative to a declin- 
ing SFR-density relation. 

Studying the origin of the SFR-density relation pro- 
vides additional insight into the star formation his- 
tories (SFHs) of galaxies across a range of environ- 
ments. For example, a declining SFR-density rela- 
tion at a fixed stellar mass can be produced by ei- 
ther (1) the proportion of quiescent and star forming 
galaxies (SFGs) changing with density, (2) the SFRs of 
SFGs decreasing at higher densities, or (3) some com- 
bination of these two scenarios. The fraction of red 
and/or blue ga laxies are often used to gauge scenario 
(1) above (e.g.. |Baldrv et al.ll2006t iCucciati et al.l 120061 : 
iGerke et al.ll2007t ICooper et al.l2007HPatel et aH20 09a) , 
however, dusty SFGs can display red op tical colors 
(|van der Wei et all 120071 iMaller et all 120091) . thus com- 
plicating such measurements. W ith extinction in galaxie s 
also dependent on environment (|Kauffmann et alll2004D . 
computing the fraction of quiescent and SFGs vs. local 
density from a single optical color is further complicated. 
Thus, selecting SFGs, both the unobscured and obscured 
variety, is essential in determining the contribution of the 
three scenarios listed above in producing the SFR-density 
relation. 

In th is work, we build on our analysis in IPatel et al.1 
(|2009bf ) in which we found declining obscured SSFR- 
density and SFR-density relations at a fixed stellar mass 
at z ~ 0.83. Here, we expand the analysis with three 
different SFR indicato rs in order to confirm these declin- 
ing relations found in IPatel et a l. (2009bJ). Our sample 
in this work also spans a l arger r edshift range (0.6 < z < 
0.9) than in IPatel et al.l (|2009bL 0.80 < z < 0.87) and 
encompasses galaxies in several groups and a cluster. As 
a consequence, we are able to determine how the SFRs of 
galaxies vary across a range of larger dark matter halos. 
Finally, one of our primary goals in this work is to deter- 
mine how the SFRs of individual galaxies change in order 



to produce the declining SFR-density relation. To this 
end, we employ a color-color diagram to distinguish qui- 
escent and SFGs, as also utilized by several other re cent 
works (e.g- lWuyts et al.ll2007t iWilliams et al.ll2009| ). 

In conducting this study, we primarily utilize stellar 
mass-limited samples as opposed to luminosity-limited 
samples. A limiting magnitude for a survey implies a 
limiting mass for old, quiescent galaxies that is generally 
higher than the limiting mass for young, SFGs at a given 
redshift (see, e.g, Appendix IA"|). As a consequence, below 
the stellar mass limit implied by passive evolution of old 
systems, no sample can be used to definitively charac- 
terize the average galaxy. This does not mean that the 
lower mass, blue SFGs found in luminosity-limited sam- 
ples cannot be used to study particular modes of star for- 
mation (SF), but merely that such samples are not neces- 
sarily representative of the full population at low mass. 
These different selections can therefore lead to widely 
varying results. For example, in studying the mass- 
limited morphology -density relation (MDR) in clusters, 
iHolden et al.l (|2007f ) found little evolution in the fraction 
of early-type galaxi es at a fixed densi ty between z ~ 1 
and z ~ (see also. IHolden et aLll2006| ). An extension of 
this ana lysis to field densities also results in a similar con- 
clusion (|van der Wei et al1l2007D . In contrast, previous 
results that utilized luminosity-limited samples found a 
significant buildup of early-type galaxies at late times, 
with the share of SOs growing at the expense o f spirals 
(e.g.. iDressler et alJll997t iPostman et all 120051 ). Other 
galaxy properties correlate with stellar mass as well, such 
as op tical col ors dBell fc de Jong 2001; Kauff mann et all 
2003), SFRs (iBrinchmann et alll2004 iSalim et al l 120071: 



Noeske et aTll2007l ). and extinction (jGarn fc Bestll2010D . 
In drawing conclusions, it is therefore important to con- 
sider the distribution of stellar masses for different sub- 
samples that are being compared. 

Our z-band spectroscopic selection (zab < 23.3 mag) 
allows us to build mass-limited sa mples at z ~ 0.9 dow n 
to M > 0.2M* (M* re 10 11 Mm. iFontana et al.ll200ll . 
probing an important stellar mass range that is respon- 
sible fo r much of the build up on the red-sequence since 
z ~ 1 (jBrown et al. 1 120071 ). This relatively low stellar 
mass limit is a key advantage of our survey over other 
large spectroscopic field and cluster surveys at these red- 
shifts. For comparison, the stellar mass limit at z ~ 0.9 
for the i?AB < 24.1 mag se lected DEEP2 sur vey is 
M re 5 X 10 10 M Q (see, e.g.. iNoeske etall 12001 . For 
the Jab < 22.5 mag selected zCOSMOS-bright survey 
the stellar mass limit is M re 7 x 10 10 M Q (see, e.g., 
iLillv et all 1 20091) . For the I < 23 mag (Vega) selected 



EDisCS survey (jWhite et al.ll2005f) the stellar mass limit 
is M re 3 x 1 10 Mm. The Jab < 2 4 mag selected VVDS- 
Deep survey (|Le Fevre et al.ll2004l ) reaches a comparable 
mass limit at z ~ 0.9 to that in this work. 

Our paper is outlined as follow. In § [2] we discuss the 
data used while in § [3] we derive key quantities such as 
stellar masses, SFRs, and rest-frame colors. We present 
our method for distinguishing quiescent and SFGs in § 3J 
Results arc presented in § [5] along with a discussion in 
§ [6] Finally, a summary is presented in § [71 

We assume a cosmology with Ho — 70 km s _1 Mpc -1 , 
Q M = 0.30, and Q A = 0.70. Stellar masses and SFRs 
are based on a Chabrier IMF (Chabrier 2003]). All mag- 



3 



nitudes are reported in the AB system unless otherwise 
stated. 



2. DATA 
2.1. Imaging 
2.1.1. Optical and Near-IR Imaging from the Ground 

The analysis in this paper is based on imaging from 
multiple telescopes and instruments that span from 
the near-UV to mi d-IR. We obta i ned a rchival Sub- 
aru Suprime-Cam (|Mivazaki et al.l I2002D imaging of 
RX J0152-13 over a wide field of size 29' x 39' (for 
a detailed discussio n of the Suprime-Cam data, see 
IKodama et aLll20"05ft . The field was observed in VRiz 
with Suprime-Cam and had a seeing FWHM~ 0'.'65 in 
all four bands. Limiting magnitudes were determined by 
placing D — 3" apertures across each image and fitting 
a Gaussian to the distribution of fluxes for the subset of 
apertures that did not fall on detected sources. The la 
width of the Gaussian fit was taken to represent the RMS 
in the background within the aperture. The resulting 5a 
depths in VRiz were 25.72, 25.65, 25.20, 24.24 AB mag 
respectively (D = 3" aperture diameter). Note that the 
z-band imaging was used for the spectroscopic selection 
and is discussed below. 

We obtained K s imaging of the central 26' x 26' using 
the W ide Field Infrared Camera (WIRC, iPersson et al.1 
2002) on the 2.5 m DuPont telescope at Las Campanas. 
The effective seeing of the mosaic was ~ 0'.'7 and the ha 
depth was 21.22 AB mag (D = 3" aperture). The K s - 
band imaging is important for constraining stellar mass 
estimates, as it samples the rest-frame J-band for most 
galaxies in our sample. When combined with the optical 
imaging, the long wavelength baseline afforded by the K s 
data also enables one to constrain the amount of dust 
extinction. 

The Magellan IMACS (|Dressler et al.l [20061 ) B-band 
imaging covers the full Suprime-Cam FOV. The seeing 
was 1'.'2 and the 5a depth was 24.55 AB mag (D = 3" 
aperture). The B-band imaging probes the rest- frame 
UV continuum and therefore provides a measure of the 
unobscured SFR. 

The BVRizK s imaging are binned to the same pixel 
scale (Of! 2), and placed on the same astrometric solution. 
The z-band image was used for source detection, and the 
"dou ble-image mode" of SExtractor (jBertin fc Arnoutsl 
1996) was used to extract fluxes in matched apertures in 
BVRizK s . 

2.1.2. Space Based Mid-IR and Optical Imaging 

We obtained Spitzer MIPS (|Rieke et al.l 12004ft 24 ^m 
imaging of the field in Cycles 2 and 5, and combined the 
data with exist ing MIPS imaging of the central regions 
of the field (see IMarcillac et al.ll2007l ). The MIPS cover- 
age in this paper with our Cycle 2 and 5 data extends 
well beyond the region centered on t he core of RX J0152 - 
13. This data was also employed in iPatel et al.l j2009b). 
The MIPS imaging covers a total o f ~ 0.1 deg 2 . We used 
MOPEX (jMakovoz fc Khadl2005l) to reduce the MIPS 
data and to create the final mosaic, which had a pixel 
scale of 1.245". The PSF FWHM was ~ 6". The APEX 
module in MOPEX was used for source detection. The 



5a detection limit was 125 /iJy for a D — 6" aperture af- 
ter applying a factor of 3.0 aperture correction to a total 
flux for point sources. MIPS sources were matched with 
the zab < 23.3 mag catalog using a l'/5 matching dis- 
tance. Our choice of matching distance reflects a trade 
off between the number of matched objects and contam- 
ination due to multiple matches. The matching distance 
is roughly equivalent to two times the RMS scatter in 
the position offset between MIPS objects and objects in 
the z-band catalog. Overall, ~ 19% of objects from the 
zab < 23.3 mag catalog with MIPS coverage are matched 
to a MIPS source. In cases where multiple objects in the 
z-band catalog could be matched with the same MIPS 
source, the closest object was assigned as the match. Of 
the objects assigned to a MIPS counterpart, ~ 4% of 
them were assigned in this way. Note that not all galax- 
ies are matched with a MIPS counterpart due to a com- 
bination of (1) galaxies with MIPS fluxes that are likely 
below the detection limit and (2) centroid shifts due to 
the confusi on limit that lead to non-matches (see, e.g., 
lHogdl200ll ). The MIP S dat a are used to derive obscured 

SFRs as discussed in § 13.71 

We also used HST ACS (|Ford et~all[200l imaging that 
was centered on RX J 0152-13. The da t a redu ction proce- 
dures are discussed in lBlakeslee et al.l (|2003l ). The imag- 
ing covers regions centered on the cluster core and out- 
skirts (R < 6')- In this paper, we use F625W and F775W 
imaging of the central regions (R < 3') and F606W and 
F814W imaging of the outskirts to make postage stamps 
(see §15). 

2.2. Spectroscopy 

We selected galaxies with z-band MAG_AUTO magni- 
tudes of zab < 23.3 mag for spectroscopy with IMACS 
on the Magellan/Baade telescope. Note that this mag- 
nitude limit is much brighter than the z-band detection 
limit. In place of the grating, we utilized a low disper- 
sion prism (LDP) designe d by S. Buries fo r use by the 
PRIMUS redshift survey (jCoil et al.l [2010h . This con- 
figuration allowed us to place ~ 3000 slits, each with a 
width of ~ 0'.'8, onto a single slitmask. The wavelength 
range of the LDP spectra is ~ 4500 A to ~ 1 /im, and 
spans roughly ~ 100 pixels. The dispersion and reso- 
lution of the LDP vary strongly with wavelength. At 
5000 A, the dispersion is ~ 20 A/pixel and the spec- 
tral resolution FWHM ~ 80 A, while at 8500 A, the 
dispersion is ~ 150 A/pixel and the spectral resolution 
FWHM - 600 A. The median total exposure time for 
each object was ~ 3 hr/pixel, roughly three times longer 
than PRIMUS. The spectra were flux calibrated using 
spectrophotometric flux standards. The data reduction 
and spectral extraction process for the LDP data are 
discussed in more detail in Patcl et al. (2011, in prepa- 
ration) . 

2.3. Data Preparation for SED Fitting 

In preparation for fitting galaxy SEDs, which include 
both broadband photometry and prism spectroscopy, we 
took steps to account for differences between the imag- 
ing PSFs. For each object, we used a D — 3" color 
aperture for the broadband photometry, which was aper- 
ture corrected to a total magnitude for all bands us- 
ing the difference between the z-band MAG_AUTO and 
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Figure 1. Difference between LDP redshifts and redshifts from 
higher resolution spectroscopy (A z = zldp — Zspec), as a frac- 
tion of (1 + Zapec) versus Zspec (top). Only the 377 galaxies with 
0.6 < z sp cc < 0.9 and zab < 23.3 mag are shown. The bottom 
panel shows a histogram of A z /(1 + z spcc ) for this sample. The 
lcr biweight scatter of A z /(1 + z S pec) is it = 1.2%. The proportion 
of catastrophic outliers is very low, with only 1.3% of the sample 
having |A Z |/(1 + z sp oc) > 5%. 

D = 3" z-band magnitude. While the VRizK s imag- 
ing had roughly similar seeing, the B-band seeing was 
slightly higher at FWHM ~ I'll. To account for the 
larger B-band PSF, we convolved the V imaging to the 
.B-band seeing of 1'.'2 and determined the difference in 
magnitude for the D = 3" aperture between the blurred 
and un-blurred U-band imaging. This difference (typ- 
ically ~ 0.04 mag) was added to the D = 3" B-band 
magnitude to "correct" it to the VRizK s seeing of ~ Of! 7. 

The LDP spectra were scaled to match the VRi pho- 
tometry. The scaling at each LDP pixel was a combi- 
nation of a constant value plus a wavelength dependent 
component, introduced to further improve upon the flux 
calibration computed from the spectrophotometric stan- 
dard. 

3. ANALYSIS 

3.1. LDP Redshifts 

LDP redshifts are determined by fitting the SEDs of 
galaxies with stellar population synthesis models. Both 
the LDP spectra and broadband photometry are used in 
the fitting. The method is briefly discussed in lPatel et al.l 



(|2009af) and will be discussed in more detail in Patel et al. 
(2011, in preparation). 

Our survey with the LDP finds galaxies spanning a 
range of redshifts. The analysis in this work focuses on 
the redshift interval 0.6 < z < 0.9. The low redshift 
boundary was determined such that the bright magni- 
tude selection criterion (zab > 18 mag) would not re- 
sult in the loss of massive galaxies. Meanwhile, the high 
redshift boundary was mostly determined by the desired 
stellar mass limit of our sample, with some considera- 
tion also for the larger redshift uncertainties at higher 
redshifts. 

To determine the precision of the LDP redshifts, zldp, 
we compare them to redshifts from hig her resolution 
spectr oscop y, z sr , Rn , using t he catalogs of iTanaka et all 
(2006) and IDemarco et ail ([20051 ). We also include 
redshifts from our own Magellan /IMACS/LDSS3 and 
Keck/DEIMOS (|Faber et al.ll2003D spectroscopy in com- 
paring zldp to z spoc in Figure [TJ For galaxies in our 
redshift range, 0.6 < z spC c < 0.9, the LDP redshifts have 
a biweight scatter of cr ZLDP _ ZBpcc /(l + z spcc ) = 1.2%. The 
scatter is approximately the same for blue and red galax- 
ies. Red galaxies have slightly lower zldp compared to 
•^spec (< 1% of (1 + z)), accounting for the non-zero value 
of the mode in Figure [TJ 

In reporting parameters from SED fitting in § 13.31 we 
used a different set of templates from those used for de- 
termining the redshift. The templates used for determin- 
ing redshifts utilized specific priors on various compo- 
nents (e.g. M/L, rest-frame colors, emission line ratios, 
etc.) in order to minimize the scatter in zldp — z spcc . 
Thus, these templates were used to optimize redshift 
measurements. In addition, new data products (BK S ) 
have been incorporated into our analysis since the red- 
shift determination process, and we have therefore em- 
ployed more flexible SFHs in fitting these data (see Fig- 
ure HJ). 

3.2. Spectroscopic Completeness 

Understanding the completeness of a spectroscopic sur- 
vey is critical in reliably computing various quantities, 
such as local galaxy densities or the typical SFR or 
fraction of objects of a certain type in a given density 
bin. We determine the spectroscopic completeness of 
our zab < 23.3 mag selected survey as a function of 
the z-band magnitude, R — z color, and position on the 
sky (see Figure [15] in the Appendix). Overall we ob- 
tained a high quality LDP redshift for 9341 galaxies with 
zab < 23.3 mag. At a given position on the sky, the 
spectroscopic completeness fraction in a magnitude-color 
bin is determined by taking the number of galaxies with 
redshift measurements, zldp, and dividing by the num- 
ber of galaxies in the same magnitude-color bin from the 
zab < 23.3 mag selection catalog. The overall spectro- 
scopic completeness of the survey is ~ 74%. The com- 
pleteness is fairly uniform for bright and faint objects, 
varying from ~ 76% at zab < 22.5 mag to ~ 71% at 
zab > 22.5 mag. Near the magnitude selection limit, 
23 < zab < 23.3 mag, the completeness is slightly lower 
(~ 60 — 70%). The completeness is also uniform in color 
at ~ 74% for R - z > 1 AB mag and R - z < 1 AB mag. 
The completeness is slightly higher in the central regions 
(R < 15': ~ 81%) compared to the outer regions of the 
LDP spectroscopy (R > 15': ~ 70%). Where appropri- 
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Figure 2. Example SED fits to LDP spectroscopy (blue) and 
BVRizK B broadband photometry (orange). The best fitting red- 
shifted full resolution BC03 r-modcl is shown in gray (NOTE: the 
actual data are fit to models that are convolved with the LDP 
resolution). The Gaussian fit to [OIL1A3727 A is indicated by the 
arrow and shown in red. An HST ACS color postage stamp (5" on 
a side) is shown inside each panel. The LDP redshift, zldPi and 
redshift from higher resolution spectroscopy, z spec , are also indi- 
cated in each panel. With the exception of the MIPS 24 fim data, 
this figure showcases all of the data sets used in this paper. The 
SF activity varies from quiescent on top to very active on bottom. 

ate, measured quantities are assigned weights based on 
these completeness maps. 



3.3. SED Fitting 

In order to measure rest-fr ame properties of ou r sam- 
ple, we fit their SEDs with iBruzual fc Charlotl (|2003l 
hereafter BC03) stellar population synthesis models (low 
resolution models), using a Chabrier IMF. The 50th per- 
centile value of zldp from the redshift likelihood function 
was used as the redshift of each galaxy and the subset of 
models from our grid (see below) with this redshift were 
used in the SED fitting. Models were fit to both the 
LDP spectroscopy and BVRizK s broadband photome- 
try. Due to the IMACS detector response and uncer- 
tainty in the flux calibration at redder wavelengths, we 
utilized the LDP wavelength range of 5000 < A < 8500 A 
in the SED fitting. 



We used a grid of r-models that span a wide 
range of galaxy SFHs. Such r-m odels SFHs are com- 
monly used in fitting SEDs ( e.g.. iPapovich et al.l [2001t 
iForster Schreiber et alJ [2001 iFranx et al.l l2008h . The 
grid consists of 4 parameters: (1) 3 metallicities: O.4Z , 
Z Q , 2.5Zq, (2) 10 values for r logarithmically spaced 
between 0.1 < t < 20 Gyr, (3) 20 ages logarithmically 
spaced between 0.3 < t < 7 Gyr (for a given redshift, ages 
are limited to the subset with values less than the age 
of the universe at that epoch), and (4) 10 values for Ay 
spaced between <Ay< 2 and assuming a lCalzetti et al.1 
(2000) extinction curve. The stellar component is com- 
bined with a Gaussian representative of [OIIJA3727 A in 
a non-negative least squares (NNLS) fit to the SED. In 
carrying out the SED fitting, several important quanti- 
ties are stored from the best-fitting minimum \ 2 model, 
including the stellar mass, SFR, and [Oil] flux from the 
Gaussian emission line component. Figure [5] shows ex- 
ample fits to SEDs, which include LDP spectroscopy and 
BVRizK s broadband photometry. 

3.4. Rest-frame Colors and Magnitudes 

We use the procedure described in Rud nick et al.l 
(2003) to apply if-corrections and determine rest-frame 
colors and magnitudes. The procedure is commonly 
used by other wo r ks to apply if -correc tions (see, e.g., 
iTavlor et al.l[200l IWflhams et al.1 12001 . To determine 
the magnitude of an object in a desired rest-frame fil- 
ter, this method uses the redshift to find the two closest 
observed filters. The observed color of the object from 
these two filters is used to search the grid of models from 
the SED fitting for the two models with the closest colors 
that straddle the observed color. We interpolate in these 
two models, using the observed color, to determine the 
magnitude of the object in the desired rest-frame filter. 
This procedure produces similar col ors and magnitudes 
at z ~ 0.8 as in IPatel et all (|2009alfbT ) . but has the ad- 
vantage of allowing us to easily expand the analysis to 
suit the larger redshift range studied in this work. The 
rest-frame filters relevant to this work are UVJ. We use 
the Buser UV, and 2MASS J filter curves supplied with 
BC03 in computing rest-frame magnitudes. 

3.5. Stellar Masses 

Stellar masses are one of the key parameters com- 
puted in the SED fitting. We primarily use a stellar 
mass limited sample to conduct our study in this pa- 
per. Given the z-band spectroscopic selection, in our 
redshift interval of interest, the limiting stellar mass is 
determined by the masses of the faint red galaxies. Al- 
though the z-band magnitude limit results in blue galax- 
ies with lower masses, their red counterparts of the same 
low mass do not make it into the sample. We illustrate 
and discuss this further in the Appendix. We find the 
limiting stellar mass at z = 0.9 is M > 1.8 x 10 10 M© 
(log M/Mq > 10.25). Of the 3326 galaxies in the red- 
shift interval 0.6 < z < 0.9 with z AB < 23.3 mag, 1174 
are above the mass limit. Of these, 977 have if s -band 
imaging. About ~ 16% of this mass-limited sample with 
if s -band imaging is located in the vicinity of the cluster 
RX J0152-13 (0.80 < z < 0.87 and R < 3 Mpc). 

3.6. Local Density 
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Figure 3. Spatial distribution of galaxies with zldp in three dif- 
ferent redshift slices, 0.6 < z < 0.7 (top), 0.7 < z < 0.8 (middle), 
and 0.8 < z < 0.9 (bottom). Galaxies above log M/Mq > 10.25 
are color-coded by local density, with divisions for the four density 
bins at 3, 13, 50 Mpc — 2 indicated by the blue, green, orange, and 
red points. Galaxies below this mass limit are shown in gray. The 
bar at the top right of each panel indicates a projected length of 
1 Mpc. An outline of the Spitzer MIPS 24 fim imaging with expo- 
sure level above ~ 1500 s is shown in black. The cluster RX J0152- 
13 is at z ~ 0.83 and is indicated by the two solid black circles that 
represent the virial radii of its two cores. The dashed circle repre- 
sents a projected clustercentric radius of 3 Mpc. 



We compute the projected local galaxy density using 
a similar but slightly updated procedure as carried out 
in lPatel et al.l (|2009bD . For a given galaxy at redshift zo, 
we determine the distance to the 5th nearest neighbor, 
g?5, where only those galaxies above the mass- limit of 
logM/M > 10.25 and with |z ncigh - z \/(l + z ) < 2% 
are included in the neighbor list. Note that this veloc- 
ity window is twice as large as the typical uncertainty 
in (1 + zldp)- In determining ds, distances to all galax- 
ies in the neighbor list are sorted from lowest to high- 
est and the corresponding weights (computed from Fig- 
ure [IS] in the Appendix) of the neighbors from the com- 
pleteness map are summed. We interpolate to find the 
distance corresponding to a sum of five in the weights 
(i.e., d§). In this way, our densities reflect values for a 
survey with 100% spectroscopic completeness. The dis- 
tance, g?5, is used to define the circular areal element 
for computing the local density. Thus, the local density 
is £ = 5/(irdf). The analysis in this paper is limited 
to the central regions of our spectroscopic coverage (i.e., 
where there is K s or MIPS imaging), therefore minimiz- 
ing edge effects. In addition, the list of neighbors for 
objects near the edge of the redshift window containing 
our sample includes objects that lie slightly beyond the 
redshift boundaries. The median density of our full mass- 
limited sample at 0.6 < z < 0.9 with i^-band imaging 
is ~ 6.5 Mpc -2 . Ignoring galaxies in the redshift inter- 
val containing the cluster, 0.80 < z < 0.87, the median 
density is ~ 4.8 Mpc -2 . 

Figure [3] shows galaxies above the mass limit in three 
redshift slices and color coded by local density. The clus- 
ter RX J0152-13 is at a redshift of z ~ 0.83 and is easily 
seen in the center of the bottom panel. Note how this red- 
shift slice, 0.8 < z < 0.9, contains a significant amount 
of structure compared to the two lower redshift slices. 
Another prominent over density is the l arge s tructure at 
z ~ 0.75. We note that iTanaka et ail (|2006f) identify it 
as a cluster but measure a velocity dispersion for this 
structure of only 210 ± 98 km s -1 from a sparse sample 
of objects. The nature of this overdensity is therefore 
unclear. 

3.7. Star Formation Rate Indicators 

We use three different measures of the SFR in this pa- 
per: (1) rest-frame 12-15 /im luminosities from Spitzer 
MIPS 24 ^im imaging, (2) the SFR of the best-fit model 
from the SED fitting, and (3) [OII]A3727 A line lumi- 
nosities. Each of these SFR indic ators has associa ted 
with it various systematics (see e.g.. lKennicuttl[l998() . It 
is therefore helpful to use multiple indicators in drawing 
conclusions about the SF activity of galaxies in different 
environments. 

3.7.1. Spitzer MIPS 24 fim 

The mid-IR emission traced by MIPS in our redshift 
interval of study (rest-frame 12 — 15 /im for our sam- 
ple) represents a combination of thermal emission from 
dust grains that reprocess UV light as well as emis- 
sion lines fro m polycyclic aromatic hydrocarbons (PAHs, 
see review bv lPuget fc Legen[T989T ). This mid-IR lumi- 
nosity correl ates with the total in frared luminosity, Lir 
(8 — 1000 /im, ICharv fc E lbaz 2001), which serves a tracer 
for the total amount of re-emitted UV- visible light. Typ- 
ical conversions between Ltr and SFR are calibrated 
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with starbursts that have ages of < 10 8 yr (K enniciittl 
1998). However, the SF timescales probed by Lir can 
vary depending on the heating source (e.g., young vs. 
old stars) and optical depth of the dust. Some works 
find that MIPS derived SFRs may represent SF over 
long timescales (~ 1 — 2 Gyr, see, e.g,. ISalim et alj|2009t 
IKelson fc HoMenl[20ll given the non-negligible contri- 
bution of dust heating associated with intermediate age 
stellar populations and TP-AGB stars. This has the ef- 
fect of overestimating the amount of SF at the epoch of 
observati on given that studies of the SF histo ry of the 
universe ([Lilly et al.1 119961: IMadau et all 119961 ) indicate 
that SFRs were higher at earlier epochs. Meanwhile, 
recent Herschel observations indicate that the observed 
MIPS 24 /im flux by itself underestimates the actual Lir 
at these redshifts when compared to Lir derived from 
the P ACS 100 /im and 160 /im bands (jRodighiero et all 
2010). The lack of consensus in the literature with re- 
gards to SFR measurements derived solely from MIPS 
24 /im data is clear. For this reason, we focus on the rel- 
ative changes in the MIPS derived SFRs across different 
environments. 

We foll ow the method of IWuvts et al.1 (|2008D and oth- 
ers (e.g., Damen et al.l 120091 ) in deriving a total infrared 
luminosity, Lir, from t he MIPS data . Thi s technique 
uses the IR templates of iDale fc Heloul ((2002) to convert 
the observed 24 fim flux into Lir. For a given redshift, 
Lir is computed for IR templates spanning a range of 
heating levels of the interstellar environment. The mean 
of all log(L IR ) is used to represent the final Lir. We note 
th at the derivation of Lir in this work differs from that 
in lPatel et all (l2009bl). To convert L m into a SFR we use 
Equation (4) from iKennicuttl (|1998D . but adjusted for a 
Chabrier IMF by multiplying by a factor of 0.56. While 
AGN also emit at these mid-IR wavelengths, they do not 
contribute significantly to the overall mid -IR luminosity 
density at these redshifts ([Bell et al.l l2005). Their contri- 
bution to the mid-IR flux for in dividual galaxies is also 
not very high (jSalim et al.ll2009l ). 

In this work, 638 galaxies at 0.6 < z < 0.9 with K s 
imaging and stellar masses M > 1.8 x 10 10 Mq have 
MIPS coverage. The depth of the MIPS imaging pre- 
cludes a galaxy-by-galaxy analysis as only ~ 28% of 
galaxies above the mass limit are detected above 75 /iJy. 
In addition, this detection rate varies for galaxies in dif- 
ferent stellar mass and local density ranges. As a con- 
sequence, it is not possible to characterize the distribu- 
tion of individual MIPS fluxes for a given sample with 
a mean or median, as many galaxies are undetected. In 
order to compute the median 24 /jm flux for a given sam- 
ple, we therefore stacked MIPS imaging for all galaxies 
in the sample, combining both detected and undetected 
sources. This approach allowed us to carry out a uniform 
analysis across various sub-samples (e.g., mass and den- 
sity) and enabled us to reach flux levels below the MIPS 
detection limit. To compute the median MIPS flux and 
its uncertainty for a given stack, we carry out a boot- 
strap analysis as follows. For a given sub-sample of N 
galaxies, we extract postage stamps from the background 
subtracted MIPS imaging, centered at the RA and Dec 
of each galaxy. We draw N random postage stamps, 
with replacement, taking into account the weights from 
the completeness map, and compute the median of the 
stack at each pixel position, thus creating a median im- 



age from the N postage stamps. Aperture photometry 
is performed on the median image using an aperture of 
diameter D — 6". The resulting flux is aperture cor- 
rected and the value stored. This process is repeated 
1000 times. The mean and standard deviation of the 
1000 measurements for the median are taken to represent 
the median flux and its uncertainty. These flux values are 
then converted into a SFR using the median redshift of 
the sample to deter mine Lir . This stac king technique 
was also employed in lPatel et al.l (l2009bl) . In thi s work, 
we build on the analysis from lPatel et al.l (|2009bj ) by ex- 
panding our sample redshift range from 0.80 < z < 0.87 
to 0.6 < z < 0.9. 

We compute the median of the stack rather than the 
mean in order to minimize the contribution from neigh- 
boring bright sources. We note that using the mean 
would have resulted in fluxes that were roughly 0.2 dex 
larger than what we found for the median. However, the 
final results would have been qualitatively the same. 

In order to test the robustness of the median flux esti- 
mates from the stacking analysis, we stacked MIPS imag- 
ing of detected objects and found the median stacked flux 
to agree with the median of the individual flux measure- 
ments to within < 10%. This was also true for detected 
MIPS objects in the highest density regions of our sam- 
ple. 

3.7.2. Star Formation Rates from SED Fit 

Each BC03 r-model that is fit to an SED has associated 
with it an instantaneous SFR. The SFRs derived from the 
SED fitting primarily leverage the rest-frame UV contin- 
uum, sampled by the observed BV filters. Stars above 
~ 5 Mq dominate the light output at these wavelengths 
and the SFR measurem ents are valid for timescales of 
10 8 yr (|Kennicuttl ll 998) . Because dust extinction is in- 
cluded in the r-models, the SED SFRs represent an ex- 
tinction corrected SFR. We caution that degeneracies 
in the parameters of the SED fitting often lead to a 
range of SFR values for a single galaxy, often depen- 
dent on the value of Ay • The typical uncertainty for Ay 
is 0.3 — 0.4 mag and was computed by refitting SEDs 
for a subset of our mass-limited sample after injecting 
Gaussian noise into the observed flux measurements and 
comparing the resulting values for Ay with the original 
ones. We use the SFR from the best fitting r-model to 
represent the SED SFR. 

3.7.3. [OII]\3727 A Emission Line Luminosities 

Nebular emission lines serve as a tracer for the SFR 
as they re-emit UV light blueward of the Lyman limit 
produced by st ars with masses o f > 10 M© and lifetimes 
of ~ 20 Myr (IKennicuttl ll998[ K As part of the SED 
fitting procedure, we include a Gaussian component for 
[OIIJA3727 A (see Figure[2|), a prominent SFR tracer that 
lies within the LDP wavelength range for our redshift in- 
terval of study. The Gaussian [OIIJA3727 A component 
is normalized such that the coefficient in the NNLS fit 
to [Oil] represents the observed line flux, Ljon] (units: 
erg s _1 cm -2 ). This line flux is converted i nto a lumi- 
nosity and then a SFR using Equation 3 from IKennicuttl 
(1998), and scaled by a factor of 0.56 for a Chabrier IMF. 

We estimate extinction corrected [Oil] luminosities 
by utilizing the best fitting Ay from the SED fit and 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 
V-J (AB mog) V-J (AB mag) 



Figure 4. Rest-frame U — V versus V — J for galaxies at 0.6 < z < 0.9. (a) UVJ diagram for a luminosity limited sample, My < —20.5 AB 
mag (Ly ~ 1.35 X 10 10 Lq). The [Williams et al. (2009) boundary separating quiescent (upper left) and SFGs (lower right) is marked 
by the gray wedge. The arrow indicates the reddening vector for AAy = 1 mag. Evolutionary tracks are shown for various BC03 solar 
metallicity SFHs. SSP (red) and CSF (blue) models are shown with time steps indicated by the diamonds at t = 0.1, 0.5, 1, 2, 3, 5 Gyr. 
After ~ 3 — 5 Gyr, the SSP falls near the quiescent clump, while the CSF lies near very blue SFGs. The colors of SFGs are extended to 
redder colors with the addition of dust. Green and cyan dashed lines represent the evolution of 1% and 5% bursts added to a t = 4 Gyr 
SSP and followed at time steps after the burst, r-ABi of *AB = 0.01, 0.1, 0.3, 0.5 Gyr (marked by squares; i-ab increases from bottom to 
top). The orange dashed line represents the 1% burst scenario but with the burst component extincted by Ay = 1 mag. The yellow star 
represents low levels of continuous SF (e.g., analogous to the Milky Way) with a CSF that has been ongoing for the last 1 Gyr combined 
with an SSP of t = 4 Gyr. The CSF component contributes 5% of the total mass. The magenta dash-dot line represents a twice solar 
metallicity r-model with r = 1 Gyr and the same time steps as the SSP and CSF marked with diamonds. In general, galaxies with any 
recent SF activity will be classified as SFGs, while those without will be classified as quiescent, (b) The subset of My < —20.5 AB mag 
galaxies with MIPS coverage. Those with MIPS detections (> 75 fJ.Jj) are indicated in red and fall primarily in the region occupied by 
SFGs. 

a ICalzetti et al.l (|2000f ) extinction law. We assume 
E(B - V)„ ns = E (B - V) stars /0A4, as indicated in 
ICalzetti et all (|20Q0[ ). We caution that these extinction 
corrections have large uncertainties given the uncertainty 
in Ay (see previous section). 

In the sections to follow, we discuss SED and [Oil] 
derived SFRs for the 977 galaxies at 0.6 < z < 0.9 with 
K s imaging and stellar masses above log M/Mq > 10.25. 
Approximately ~ 82% of galaxies in this sample have a 
non-zero [Oil] flux from the NNLS fit. Note that this 
sample is larger than the MIPS sample. 



4. GALAXY CLASSIFICATION WITH A COLOR-COLOR 
DIAGRAM: DISTINGUISHING QUIESCENT AND STAR 
FORMING GALAXIES 

One of the main science goals of this paper is to de- 
termine how the SFRs of galaxies vary across a range 
of environments. In particular, we would like to know 
how SFGs c ontribu te to the SSFR-density relation. In 
IPatel et al.l (|2009bl ) we hypothesized that the change in 
the median SFR with density was driven simply by a 
changing proportion of quiescent and SFGs, as deter- 
mined from the fraction of red galaxies, leaving a negli- 
gible contribution from any change in the SFRs of SFGs 



with density. Many other works also routinely use the 
red-sequence, as defined with a rest-frame color and mag- 
nitude (or stellar mass) , to represent quiescent galax- 
ies (e .g., iBell et al.|[200l iFaber et alj|2?j07t iBrown et al.l 
I2007T ) . However, dust can lead to a significant fraction 
of SFGs occupying the red-sequence (Ivan der Wei et al.l 
[20071 IMaller et all 12001 IWhitakcr et al.l 120101) . Below, 
we utilize a color-color selection that allows us to iden- 
tify these reddened SFGs, in addition to blue SFGs, and 
distinguish both types from quiescent galaxies. 

4.1. The UVJ Diagram 

Figure 2] shows rest-frame U — V versus V — J (here- 
after referred to as the UVJ diagram). Such color-color 
diagrams have been utilized by several ot her works to 
distin g uish SFGs from quies c ent ga l axies dWuvts et all 
1 20071 IWilliams et all [2009T [2010: IWolf et all 120091: 
iBalogh et al.H2009t IWhitaker et al.ll2010f ). While U - V 
primarily identifies age, V — J is sensitive to age and red- 
dening owing to the long baseline in wavelength. One key 
feature of the UVJ diagram is the red "quiescent clump" 
(JJ - v ~ 2 AB mag and V - J ~ 1.3 AB mag), the sub- 
sample of red-sequence galaxies that lack SF, as shown 
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later. Below the quiescent clump lies a se quence of SFGs. 
We us e the color boundaries derived by William s et al.l 
( 2009) to separate quiescent galaxies from SFGs. Model 
BC03 evolutionary tracks are shown for an SSP (red line) 
and "constant" SFR (CSF, blue line) model (r-model 
with r = 20 Gyr) for solar metallicity. The SSP model 
track passes near the quiescent clump after t = 3 Gyr, 
while the CSF model track terminates at the blue end of 
where SFGs lie. The reddening vector in Figure|4ji shows 
that SFGs are extended to redder colors with the addi- 
tion of dust extinction. Note that some of the reddest 
SFGs (V — J > 1.6 AB mag) would have been classified 
as being "red-and-dead" based on U — V colors alone. 
Instead, these galaxies are distinct from the quiescent 
clump and are predicted to be SFGs with high values 
of Ay. Bicolor diagrams, such as the UVJ ones shown 
here, can therefore be used to identify SFGs, including 
those that are reddened. 

4.2. A Diversity of Star Formation Histories for 
UVJ -selected SFGs 

While the SSP and CSF model tracks in FigurelUz pro- 
vide a glimpse of how simple galaxy SFHs appear in the 
UVJ diagram, the actual SFHs of galaxies are likely to be 
more complex. In FigureHJi we also show how a variety of 
SFHs would appear in the UVJ diagram and emphasize 
the effectiveness of the Williams et al.l (|2009f) boundary 
in distinguishing quiescent and SFGs. These examples 
also demonstrate that SFGs need not have formed the 
bulk of their mass recently in order to be classified as an 
SFG. Even the slightest amount of recent SF can move 
formerly quiescent galaxies into the SFG region of the 
UVJ diagram. 

SFHs that proceed with several ongoing bursts have 
been suggested as a mode for galaxies to form a large 
fraction o f their stellar mass, especially at th ese higher 
redshifts (|Bell et al.ll2005t iDressler et al.ll2009h . The tra- 
jectories of various bursting SFHs are shown in FigureHJi 
as dashed lines. The green and cyan colored tracks show 
the evolution of a 1% and 5% burst (by mass) added 
to a 4 Gyr solar metallicity SSP at time steps after the 
burst, £ab, of 0.01, 0.1, 0.3, and 0.5 Gyr (time steps 
shown as squares). Immediately after the burst (i.e., 
*ab = 0.01 Gyr), the U — V colors become bluer then 
most SFGs in our sample. After £ab ~ 0.1 — 0.3 Gyr, 
the colors once again return to resemble a galaxy re- 
siding in the quiescent clump. This transition proceeds 
more quickly for the smaller 1% burst. These models also 
demonstrate that unreddened bursts on top of old stel- 
lar populations are likely not capable of producing SFGs 
with the reddest colors in the UVJ diagram. Those red 
SFGs are most likely to be obscured at some level. 

The orange track in Figure |4] represents a 1% burst 
added to an SSP but with the burst component extincted 
by Ay = 1 mag. Note the displacement from the unob- 
scured 1% burst (green track). These dusty bursts spend 
even less time (< 0.1 Gyr) classified as SFGs. Depend- 
ing on the level of extinction and the distribution of the 
dust, reddened bursts can populate the reddest colors for 
SFGs. 

A galaxy with low levels of continuous SF, qualitatively 
resembling the state of the Milky Way, is represented 
by the single yellow star at U — V ~ 1.2 AB mag and 
V — J ~ 0.9 AB mag. This SFH is a combination of a 



CSF that has been ongoing for the last 1 Gyr and an SSP 
of age 4 Gyr. The CSF component has produced 5% of 
the total mass. Such low levels of continuous SF would 
place a galaxy in the region occupied by SFGs. 

SFHs with SFRs that slowly decline with time can lead 
to galaxies being classified as SFGs despite very low SS- 
FRs, especially for super-solar metallicity stellar popu- 
lations. The magenta dash-dot track represents a twice 
solar metallicity r-model with r = 1 Gyr. Time steps 
similar to that of the SSP and CSF are shown for this 
model. The use of a super solar metallicity model gener- 
ally results in redder V — J colors at fixed U — V . After 
3 — 5 e-folding times in the SFR (i.e., between the last 
two time steps), this r-model would remain classified as 
a SFG with a SSFR between - 0.1 and 0.02 Gyr" 1 . We 
note that for the corresponding solar metallicity SFH 
(not shown), the model track would cross the boundary 
between SFG and quiescent after ~ 3 Gyr. Thus, some of 
the galaxies in the "green valley" (1.2>f/-F>1.6AB 
mag) are likely to be comprised of galaxies with some- 
what older stellar populations with some ongoing SF, 
rather than SFGs with more vigorous SF and modest 
amounts of dust (e.g., Ay ~ 1 mag). 

To summarize, based on BC03 stellar population syn- 
thesis models, galaxies classified as star-forming in the 
UVJ diagram can be comprised of (1) galaxies with rel- 
atively constant SFRs since the onset of SF, (2) formerly 
quiescent galaxies that have undergone a very recent 
burst, (3) galaxies with some residual, low level, long- 
term SF (e.g., 5% of mass added over the last 1 Gyr 
with roughly constant SFR), and (4) super-solar metal- 
licity galaxies with low levels of SF. The addition of dust 
to any of these scenarios moves SFGs towards redder 
U — V and V — J colors, and generally keeps them in 
the region occupied by SFGs. Thus, galaxies with virtu- 
ally any level of recent SF activity would be classified as 
SFGs, while thos e without would be cl assified as quies- 
cent based on the iWilliams et al.l (|2009l ) boundary. 



4.3. MIPS Detections Trace UVJ -Selected SFGs 

The models discussed above provide a theo retically 
motivated backing for the IWilliams et al.l (|2009t ) bound- 
ary separating quiescent and SFGs. Here we perform an 
empirical check on this boundary. Figure HJ> shows galax- 
ies with MIPS 24 fjLm. coverage and indicates the subset 
with 24 fim detections (> 75 /xJy). These detections 
primarily trace the region occupied by SFGs. Note the 
detections at red colors (U — V > 1.6 AB mag) where 
galaxies would have been considered members of the red- 
sequence based solely on their U — V color. Most of these 
MIPS detections are in the region of the UVJ diagram 
consistent with SFGs with extinction. A stacking anal- 
ysis of the SFGs in this region (U — V > 1.6 AB mag) 
without MIPS detections suggests that they are slightly 
below the detection limit with a median 24 /im flux of 
~ 41 fiJy, and therefore have modest obscured SF activ- 
ity as well (~ 3 M Q yr _1 at z = 0.75). In contrast, the 
median stacked MIPS 24 fim flux for quiescent galaxies at 
U - V > 1.6 AB mag is a factor of - 3 lower. The MIPS 
detections therefore lend further support to our ability 
to isolate SFGs in the UVJ diagram, including both blue 
(i.e., unattenuated) and red (i.e., dust reddened) SFGs. 
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Figure 5. UVJ diagram color-coded by the median SED SSFR 
in color-color bins of size 0.1 mag. (a) Galaxies at 0.6 < z < 
0.9 with rest-frame My < —20.5 AB mag. (b) Galaxies above a 
stellar mass of log M/Mq > 10.25. The colors that correspond 
to different log(SSFR) are shown in the bottom right. Bins with 
fewer than three obj ects are not included in the analysis. The 
Williams ct al. (2009) boundary separating quiescent galaxies from 
SFGs is shown in black. Note how selecting by mass minimizes the 
bias against faint red galaxies and therefore results in fewer blue 
SFGs compared to the luminosity selection in (a). 

4.4. The Distribution of Star Formation in the UVJ 
Diagram 

The models tracks shown in Figure |4] provide clues on 
the recent SFHs of galaxies that occupy different regions 
of the UVJ diagram. However, it is helpful to know how 
the SFRs or SSFRs of galaxies in our sample are dis- 
tributed in the diagram, thus providing a check on our 
division between quiescent and SFGs. Figure [5] shows 
the UVJ diagram divided into color-color bins of size 
0.1 mag, with each bin color-coded by the median SSFR 
determined from the SED fit. The top panel shows a 
luminosity-selected sample of galaxies at 0.6 < z < 0.9 
with My < —20.5 AB mag, while the bottom panel 
shows a mass-selected sample with log M/Mq > 10.25. 



Note that selecting galaxies above a stellar mass limit of 
log M/Mq > 10.25 minimizes the bias against faint red 
galaxies and therefore results in fewer blue SFGs com- 
pared to a luminosity selection. Lines of constant SSFR 
span from blue U — V and V — J colors to redder ones 
for SFGs. Thus, SFGs with red U — V and V — J col- 
ors indicative of high levels of reddening, perhaps ow- 
ing to their inclination, can have similar SED based SS- 
FRs to those with bluer colors. At a fixed V — J color, 
SFGs with lower SED based SSFRs have redder U - V 
colors. The recent SFH of a galaxy therefore also de- 
termines where SFGs lie in the UVJ diagram, in addi- 
tion to the level of ex tinction. Figure [S] shows that the 
iWilliams et al.l (|2009f ) boundary between quiescent and 
SFGs is well justified. Note too that although quiescent 
galaxies with bluer U — V colors (U — V < 1.7 AB mag) 
have somewhat higher SSFRs than redder ones, there are 
very few galaxies in that region of UVJ color-space (see 
Figure SJi). 

Finally, it is worth noting that while there appears 
to be a strong gradient a t a fixed V — J color in SED 
based SSFRs for SFGs, IWilliams et all (pOOl find a 
much smaller r ange of SSFRs for MIPS based SFRs. 
IWilliams et al.l (j2010l ) however also find a stronger gra- 
dient for SED based SSFRs. These differences between 
SFR indicators are important to consider when interpret- 
ing the results below, especially for the SFRs of SFGs 
(see 



4.5. A Morphological Overview of the UVJ Diagram 

Figure HI populates the UVJ diagram with HST ACS 
color postage stamps for a luminosity-limited sample. 
These high resolution images provide a sense for how disk 
or bulge dominated systems, and their orientations, help 
to explain what is seen in different regions of the UVJ 
diagram. The quiescent region is dominated by galax- 
ies with strong bulge components while the star-forming 
region is comprised of disks. Blue SFGs appear to be 
face-on disks and the bluest of them (U — V ~ 0.7 AB 
mag) look like bursting systems - consistent with the 
location in UVJ color space predicted by the bursting 
tracks discussed in Figure 0Ji. Meanwhile, the reddest 
SFGs appear to be edge-on disks, consistent with the 
red colors arising from extinction. Interestingly, several 
SFGs with intermediate colors (U — V ~ 1.6 AB mag and 
V — J ~ 1.4 AB mag) appear to have clear disks but also 
host a strong bulge component. The colors of these sys- 
tems likely arise from older stellar populations, perhaps 
in a transition phase, rather than dust. A more detailed 
discussion of the morphological structure of galaxies in 
various environments and in the UVJ diagram will be 
presented in a future paper (Patel et al. 2011, in prepa- 
ration). 

4.6. The UVJ Diagram in Different Environments 

Our analysis in § [5] utilizes the UVJ diagram in dif- 
ferent density bins in order to examine relative propor- 
tions of quiescent and SFGs as well as the SFRs of SFGs. 
Figure [7] shows the UVJ diagram in four density bins 
for galaxies at 0.6 < z < 0.9 above the mass limit of 
log M/Mq > 10.25. Note the presence of many galax- 
ies with red U — V colors (U — V > 1.7 AB mag) that 
are classified as SFGs. For galaxies above the mass limit 
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Fig ure 6. UVJ di a gram represented by HST ACS postage stamps (3" on a side) of galaxies at 0.6 < z < 0.9 with My < —20.5 AB mag. 
The IWilliams eiT al, (2009) boundary separating quiescent (top left) and SFGs (bottom right) is marked in white. Quiescent galaxies are 
dominated by bulges. The bluest SFGs (V — J < 0.8 AB mag) are primarily disks, while the reddest SFGs (V — J > 1.6 AB mag) appear 
to be inclined disks. 



that would be classified as red-sequence members based 
solely on their U — V colors (based o n a red-blue separa - 
tion method similar to that applied in lPatel et all 2009a), 
roughly ~ 32 % are c lassified here as UV J-selected SFGs. 
iBalogh et al.l (|2009f ) find a similar value at somewhat 
lower redshifts and for a slightly lower mass limit. Thus, 
a considerable portion of optically red galaxies in our 
redshift interval are dusty and star-forming. 

Figure [7] also indicates the median and standard de- 
viation of the logarithm of the stellar mass in different 
density bins for our galaxy sample above the mass limit 
of log M/Mq > 10.25 (black), as well as for the sam- 
ple of SFGs (blue). The median stellar mass varies by 
less than < 0.15 dex between different density bins for 



both samples. The scatter in stellar mass between den- 
sity bins is also comparable. Given that the distribu- 
tion of stellar masses in different density bins are roughly 
similar, mass-dependent correlations with various galaxy 
properties cannot be solely resp onsible for any dens ity- 
dependent correlations (see also. lCooper et alTl2010T ) . 

5. A DECLINING SSFR-DENSITY RELATION AT 
0.6 < z < 0.9 

In this section, we study the SSFR-density relation 
for galaxies at 0.6 < z < 0.9 with stellar masses above 
log M/Mq > 10.25. We utilize the three different SFR 
indicators discussed in § 13.71 (1) rest-frame 12-15 /im 
luminosities from Spitzer MIPS 24 /im imaging, (2) SED 
fitting, and (3) extinction-corrected [OII]A3727 A emis- 
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Figur e 7. UVJ dia gram for galaxies at 0.6 < z < 0.9 with mass logAf/A/0 > 10.25 in four density bins. The red line denotes the 
Williams ct al. (2009) boundary separating quiescent and SFGs. Note that selecting galaxies above a mass limit of log M/Mq > 10.25 
results in fewer blue SFGs compared to the luminosity-limited selection in Figure [4] The median and standard deviation of log M/M@ for 
the full mass-limited sample (black) and of the SFG sample (blue) are labeled in the bottom right for each density bin. The mass function 
is fairly uniform in different environments for both samples, with slightly more massive galaxies at the highest densities. Mass dependent 
correlations for galaxy properties in different density bins are therefore minimized by using a stellar mass limited sample. 
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Figure 8. Median SSFR versus local density for galaxies with 
mass log M/Mq > 10.25 at 0.6 < z < 0.9. Gray triangles on 
the x-axis denote boundaries of density bins. SFRs are derived 
from three different indicators: (1) rest-frame 12-15 fim luminosi- 
ties from Spitzer MIPS 24 fim imaging (red circles), (2) SED fitting 
(green squares), and (3) extinction corrected [OII]A3727 A emis- 
sion line luminosities (blue diamonds). The colored lines indicate 
the best fitting log(SSFR)-log(density) relations with correspond- 
ing parameters given in Table [T] For all three SFR indicators the 
SSFR-density relation decreases, from the low-density field to the 
cores of groups and a cluster. 

sion line luminosities. Each of these SFR indicators has 
associated with them various systematic uncertainties. 
For the purposes of this paper, we are interested in the 
relative change in the SFR within a single indicator, 
across different density bins. We note that the MIPS 
derived SFRs are limited to the sample with MIPS cov- 
erage (see Figure [3]), while the SED and [Oil] SFRs are 
reported for a larger sample. 

Figure [5] shows the median SSFR-density relation for 
each of the three SFR indicators. Uncertainties were 
computed for the SED and [Oil] SSFRs by bootstrapping 
the sample within each density bin. A bootstrap tech- 
nique was also used to compute the uncertainties for the 
MIPS SSFRs (see § [3771). The position of the data points 
along the x-axis reflects the median value for the local 



density within the given sub-sample. At 0.6 < z < 0.9, 
the SSFR from all three SFR indicators declines across 
the entire range of densities, from the low-density field 
to the cores of groups and a cluster. Thus, the over- 
all SF activity in the cores of groups and a cluster at 
this epoch is lower than in the field. Straight line fits to 
the MIPS, SED, and [Oil] log(SSFR)-log(density) rela- 
tions indicate non-zero, negative slopes at significances 
of ~ 5, 6, and 7cr. The best fit parameters are given 
in Table [TJ The SED and [Oil] based SSFRs decline 
by roughly an order of magnitude across the full density 
range, while the MIPS based SSFRs decline by a factor 
of ~ 4. We note that the raw [Oil] SSFRs that are uncor- 
rected for extinction (not shown), which depend on the 
SED fits only in establishing a continuum level above 
which the line flux is integrated, also show a declining 
SSFR with increasing density, although with a shallower 
slope (A(log SSFR) /A (log £) w -0.35±0.05). This sug- 
gests relatively higher levels of extinction at lower densi- 
ties (see jj EH) . 

A large proportion of our galaxy sample in the highest 
density environments reside in the vicinity of the clus- 
ter. Some of these galaxies have likely passed through 
the cluster core at so me point in time, res ulting in their 
lower SFRs (see e.g.. IBalogh et all 120001) . It is there- 
fore reasonable to ask if the cluster is solely responsi- 
ble for driving the declining SSFR-density relation that 
we found in Figure [5] Figure [§] shows the mass-limited 
SSFR-density relation for all three SFR indicators af- 
ter removing galaxies near the cluster with redshifts in 
the range 0.80 < z < 0.87. Removing galaxies in this 
redshift range primarily impacts the two highest density 
bins (note the larger error bars). The SSFR-density rela- 
tion continues to decline at higher densities for all three 
SFR indicators, although the decline is not as steep for 
the MIPS based SFRs. Straight line fits to the MIPS, 
SED, and [Oil] log(SSFR)-log(density) relations yield 
non-zero, negative slopes at the ~ 2.5, 5, and 5a levels 
respectively (see Table Q] for best fit parameters). Within 
the uncertainties, the SSFRs of galaxies far from the clus- 
ter in the highest density bin are essentially the same as 
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Table 1 

Parameters for best-fit log — log lines for SSFR-density 
relation 



SFR Indicator Intercept Slope 

(log SSFR) (A(log SSFR)/A(logE)) 



All galaxies (Figure [HJ 



MIPS -9.61 ±0.082 -0.366 ±0.072 

SED -9.95 ±0.099 -0.755 ±0.12 

Oil -9.82 ± 0.078 -0.530 ± 0.073 




Figure 9. Median SSFR versus local density for galaxies with 
mass log M/Mq > 10.25 at 0.6 < z < 0.9, but only for galaxies 
far from the cluster RX J0152-13 (outside of 0.80 < z < 0.87). 
The SSFR-density relation continues to decline, indicating that 
galaxies near the cluster environment are not solely responsible for 
the trend. Galaxies in groups at these redshifts also have lower 
SSFRs compared to the field. 

when including the cluster. We note that we arrive at 
the same conclusion when also excluding the structure 
at z ~ 0.75. This demonstrates that our result of a de- 
clining SSFR-density relation in Figure [8] is general to 
a range of environments and not solely driven by galax- 
ies near a rich cluster. Galaxies in groups at these red- 
shifts therefore also have lower median SSFRs compared 
to galaxies in the field. In the sections that follow, we 
study the full mass-limited sample at 0.6 < z < 0.9 but 
emphasize that the results are qualitatively similar when 
using a sample that ignores the cluster redshift slice. 

5.1. Does Environment Matter or is Stellar Mass 
Driving the Declining SSFR-Density Relation? 

Higher density environmen ts are found to host 
more massive ga laxies (e.g., Kauffma nn et al.l l2004t 
iBaldrv et all [20061 ) . Because SFRs and SSFRs are cor- 
related with stellar mass, it is important to account for 
differences in the distribution of stellar masses between 
environments in order to determine whether the envi- 
ronment plays an independent role from stellar mass in 
impacting galaxy SFRs. Utilizing a simple mass cut, as 
shown in Figures [8] and |H to first order does a good job 
of minimizing the differences in the stellar mass distribu- 
tions between environments (see numbers in Figure [7]). 
Here, we go a step further and study the SSFR-density 
relation at a fixed mass by dividing our sample into three 
mass bins. 

Figure [TU] shows the median SSFR versus local den- 
sity for all galaxies at 0.6 < z < 0.9 in three mass 
bins for SFRs derived from SED fitting and extinction 
corrected [OII]A3727 A. We note that a si milar figure 
for the MIPS based SFRs was presented in iPatel et all 
(2009b|) and arrived at qualitatively similar conclusions 
as presented for the other two SFR indicators below. The 
mass bins in Figure [10] are the same for each SFR indi- 
cator and were partitioned such that for the full sample 
with K s imaging, the number of galaxies in each bin 
was roughly the same. The low mass bin ranges from 
10.25 < log M/Mq < 10.50, the intermediate mass bin 



Cluster removed (Figure [9j 



MIPS -9.73 ±0.096 -0.246 ±0.10 

SED -9.98 ± 0.12 -0.653 ±0.12 

Oil -9.82 ±0.094 -0.534 ±0.098 



UV J-selected star forming galaxies (Figure 1 1 2 1 > 

MIPS -9.46 ±0.086 -0.0787 ± 0.090 

SED -9.61 ±0.11 -0.472 ±0.11 

Oil -9.42 ±0.099 -0.438 ±0.10 



from 10.50 < log M/Mq < 10.80, and the high mass bin 
from log M/Mq > 10.80. The widths of the mass bins, 
especially the two lowest, are fairly narrow, allowing one 
to investigate the SSFR-density relation at a fixed mass. 
For both SFR indicators and for all mass bins, the SSFR- 
density relation declines, from the low density field to 
the cores of groups and a cluster. Table [2] indicates the 
best fit parameters to the log(SSFR)-log(density) rela- 
tions for each mass bin. The relations for all three mass 
bins have negative slopes that are significant at > 3a. 
Any variation in the stellar mass distribution within a 
mass bin across different environments would not be suf- 
ficient to explain the order of magnitude decline in SSFR 
found in Figure [TO] This confirms that the local environ- 
ment impacts the median SSFRs and SFRs of galaxies at 
0.6 < z < 0.9 when controlling for stellar mass. 

Figure [10] also shows that for the two highest mass 
bins, the spread in SSFR at a fixed density is fairly 
small, implying that the SSFR-density relations are com- 
parable. Galaxies in these two mass bins lie above 
log M/Mq > 1 0.5, which is similar to t he characteristic 
mass found by [Kauffmann ct al. (|2003T > at z ~ above 
which galaxy properties appear similar. Interestingly, 
the SSFRs for the highest mass galaxies (log M/Mq > 
10.8) at 0.6 < z < 0.9 are also sensitive to their local 
environment, indicating that stellar mass is not solely 
responsible for driving their SF properties. 

Finally, we note that because the SSFR-density rela- 
tions are shown for narrow mass bins, the SFR-density 
relations for a given mass bin have slopes that are simi- 
lar (in log-log space) to the corresponding SSFR-density 
relations. 

5.2. Quiescent Fraction at 0.6 < z < 0.9 as a Function 
of Density 

Here, we investigate how a changing proportion of qui- 
escent and SFGs with density contribute towards the de- 
clining SSFR-density relation found in § [5] Figure QT] 
shows the fraction of UV J-selected quiescent galaxies at 
0.6 < z < 0.9 as a function of local density using the 
mass-limited sample. As one might expect from the in- 
creasing red galaxy fraction (|Patel et al.ll2009"al) for this 
sample, galaxies at higher densities are more likely to be 
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Figure 10. Median SSFR versus local density for galaxies at 
0.6 < z < 0.9 in three mass bins (low - teal circles, intermedi- 
ate - orange diamonds, high - purple squares) containing roughly 
equal numbers of galaxies, (a) SFRs derived from SED fitting, (b) 
SFRs derived from extinction corrected [OII]A3727 A. The colored 
lines indicate the best fitting log(SSFR)-log(density) relations with 
corresponding parameters given in Tablel2l At a fixed stellar mass, 
the SSFR declines from low to high density for both SFR indicators 
by roughly an order of magnitude. 

in the quiescent clump. The fraction of galaxies in the 
quiescent clump increases from 32 ± 3% at low density 
to 79 ± 4% at high density. The quiescent fraction re- 
mains fairly constant across the two lowest density bins, 
which are representative of the held. The fraction begins 
to increase substantially between the middle two density 
bins. Figure [3] shows that these two density bins rep- 
resent the transition between what are likely filaments, 
small groups, and the outskirts of larger groups (green) 
to the centers of large groups and the outskirts of clusters 
(orange) . Several other works at these redshifts also find 
evidence for a lower fraction of actively star forming sys- 
tems in higher density e nvironments, ba sed on Ha emis- 
sion ([Sobral et al.ll20lil ). [O il] emission (iPoggianti et al.l 
2008j) or mid-IR emission (jFinn et al.l I2010D (or both 
emission li nes and mid-IR flux at somewhat lo wer red- 
shifts, e.g. JTran et "all [20091: iBalogh et al.ll2009l) . 

Finally, we note that whe n using a similar mass cut as 
in Ivan der Wei et~a l. (2007), the quiescent fraction fol- 
lows the early- type galaxy fraction vs. density from that 
work remarkably well. 

5.3. The Role of Star-Forming Galaxies in the 
SSFR-Density Relation 

The analysis in § [5] for all galaxies above the mass 
limit of log M/Mq > 10.25, establishes that the me- 
dian SSFR and SFR decline in higher density environ- 
ments at 0.6 < z < 0.9, confirming the declining trend 



Table 2 

Parameters for best-fit log — log lin es f or SSFR-density relations in 
Figure [10] 



Mass range 


Intercept 
(log SSFR) 


Slope 

(A(log SSFR)/A(log £)) 




SED 


10.25 < log M/Mq < 10.50 
10.50 < log M/M Q < 10.80 
10.80 < log M/Mq 


-9.55 ±0.19 
-10.1 ±0.20 
-10.2 ±0.12 


-0.810 ±0.19 
-0.721 ±0.16 
-0.674 ±0.13 


Oil 


10.25 < log M/Mq < 10.50 
10.50 < log M/Mq < 10.80 
10.80 < log M/Mq 


-9.36 ±0.14 
-9.87 ±0.14 
-10.1 ±0.11 


-0.660 ±0.12 
-0.499 ±0.13 
-0.375 ±0.11 
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Figure 11. Fraction of UV J-selected quiescent galaxies versus 
density for galaxies with mass log M/Mq > 10.25 at 0.6 < 2 < 0.9. 
The fraction of galaxies in the quiescent clump more than doubles 
from 32% at the lowest densities to 79% in the cores of groups and 
a cluster. 

found in [Fatel et al.1 (|2009H ). In the previous section, 
we found that the changing proportion of quiescent and 
SFGs contributes towards the declining SF activity at 
higher densities. Here, we investigate whether a change 
in the SFRs of UV J-selected SFGs contributes towards 
the declining SSFR-density relation. 

Utilizing the UV J diagram (Figure [7]) , we now select 
SFGs. Figure H2] shows the SSFR-density relation, de- 
rived from all three SFR indicators, for SFGs above the 
mass limit at 0.6 < z < 0.9 with K s imaging. The 
SSFRs of SFGs decline at higher densities for the SED 
and [Oil] derived SFRs by factors of ~ 5 — 6, while the 
MIPS derived SSFRs decline only by a factor of ~ 1.3 
across the full range of densities. Straight line fits to 
the MIPS, SED, and [Oil] log(SSFR)-log(density) rela- 
tions indicate non-zero, negative slopes at significances 
of ~ 0.9, 4, and 4c (see Table [1] for best fit parameters). 
We note that if galaxies at 0.80 < z < 0.87 are excluded 
(i.e., those near the cluster), the MIPS SSFRs for SFGs 
decline by a factor of ~ 1.5 from the lowest to highest 
densities but the negative slope from a straight line fit 
(A(log SSFR) /A (log S) w -0.10) is only significant at 
the ~ l.Ocr level. The decline in the SSFR-density rela- 
tion, especially for the SED and [Oil] derived relations, 
is therefore driven in part by declining SFRs of SFGs, 
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Figure 12. Median SSFR versus local density for SFGs with mass 
log M/Mq > 10.25 at 0.6 < z < 0.9. All three SFR indicators find 
the SSFRs of SFGs to decrease at higher densities. The decline in 
the SSFR-density relation for all galaxies (Figure [8j is therefore 
caused by declining SFRs for SFGs as well as a change in the pro- 
portion of quiescent and SFGs (Figure 1 11 D . Note that the decline 
in SSFR with density for the MIPS derived SFRs is not as strong 
as the other two SFR indicators. 

in addition to the increasing fraction of quiescent galax- 
ies at higher densities. We note that at a given density 
the median SSFRs for quiescent galaxies are typically 
more than an order of magnitude lower than the value 
for SFGs for a particular SFR indicator, and generally 
lack any significant dependence on density. 

5.4. A Parallel Decline with Density in Star- Formation 
and Dust Content 

SF and dust production are connected since most dust 
is injected into the ISM through AGB stars or super- 
nova e, both the by -products of recent episodes of SF (see, 
e.g., IMathisI [19901) . Subsequent episodes of SF are also 
dependent on the presence of dust in allowing gas to cool 
and form stars. With the SFRs of SFGs decreasing with 
increasing density, we should expect to see less dust as- 
sociated with SFGs at higher densities. 

Figure [T3b shows the extinction corrected [Oil] SFR- 
density relation for the same sample of SFGs used in Fig- 
ure [T2] A straight line fit to the log(SFR)-log(density) 
relation gives a slope of A (log SFR)/A(logS) « -0.33 
and is significant at ~ 3.6cr. The decline in SFR from 
low to high density is a factor of ~ 3, somewhat smaller 
than the factor of ~ 5 decline in SSFR (Figure [T2]) . How- 
ever, within the uncertainties on the best-fit slopes, these 
numbers are compatible. Note also that the median mass 
is slightly higher at higher densities (see blue numbers 
in Figure [7]), leading to a steeper slope for the SSFR- 
density relation. The [Oil] derived SFRs at high den- 
sities co uld also be overestimated due to AGN ac tivity 
(see e.g.. lLemaux et al] [2010; Koccvski ct al. 2010b), fur- 
ther strengthening the declining SSFR-density and SFR- 
density trends found here. 

Meanwhile, Figure 113b shows the median value of the 
best-fitting Ay from the SED fits versus density for the 
same SFGs. The lowest density environments at 0.6 < 
z < 0.9 are associated with the highest dust content, with 
the median Ay decreasing by ~ 0.5 mag from low to high 
density. A straight line fit to the Ay-log(density) relation 
indicates a slope of A(^4y)/A(logS) ~ —0.32 mag dex -1 
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Figure 13. (a) Median extinction corrected [Oil] SFR versus local 
density for SFGs with mass log M/Mq > 10.25 at 0.6 < z < 0.9. 
(b) Median Av from SED fit versus local density for the same 
sample of SFGs. Both the attenuating dust content and SFRs of 
SFGs decline at higher densities, as one would expect if the cold 
gas supply is diminished in such environments. 

and is significant at the ~ 3.0c level. We note that the 
median Ay of quiescent galaxies remains constant across 
all densities, suggesting that the relative differences in 
Ay for SFGs are credible. The parallel decline in SFR 
and Ay with density for SFGs suggests that the cessa- 
tion of SF in higher density environments is linked with 
the level of attenuating dust. The lower dust content at 
higher density could be a consequence of a physical mech- 
anism removing the dust, or a lack of replenishment due 
to lower levels of SF. 

6. DISCUSSION 

6.1. No Reversal in the Mass-limited SSFR-Density and 
SFR-Density Relations at z < 1 

At 0.6 < z < 0.9, we find that at a fixed stellar mass, 
the median SSFR and SFR declines from low to high den- 
sity environments for three differen t SFR indicators, a s 
also found with MIPS based SFRs in lPatel et all (|2009bl) . 
Another new key result is that we find this declining 
SFR-density relation persists even when excluding galax- 
ies in the vicinity of the cluster (Figure|H]), indicating that 
the declining SFR trend with densit y applies mo r e gener - 
ally to all environments. As note d in|Patel et"aLl (120 09^. 
our result differs from that of lElbaz et al.l (|2007t ) and 
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I Cooper etT al, (2008), both of which find a reversal in the 
SFR-density relation at z ~ 1 for luminosity-limited sam- 
ples. The reversal found by these other works does not 
extend to the highest densities probed by their surveys. 
Instead, their SFR-density relations turnover and de- 
crease at intermediate densities and above. We note how- 
ever that we do not find any intermediate density regime 
in which the median SFR is significantly elevated rela- 
tive to the field. Our results indicate that for galaxies at 
0.6 < z < 0.9 with stellar masses M > 1.8 x 10 10 M©, SF 
activity in groups and clusters is lower than i n the field, 
simil ar to what is found at lower redshift s (jWolf et al.l 
I2009D and at z ~ (|Kauffmann et alj|2004h . Other work 
at z ~ 1 also find a declin ing SFR or SSFR-density rela- 
tion (|Scoville et al.ll200l . 

The size of the density bins used in this study are fairly 
broad (~ 0.5 dex). Is it possible that at the intermediate 
densities where a reversal has been suggested, our den- 
sity bins lack the resolution requir ed to sample any ele - 
vated SF activity? The reversal in iCooper et al.l (2008) 
occurs at overdensities relative to the median density of 
~ 0.06 to ~ 5. This translates into a density range that 
rough ly covers our three lowest density bins. The rever- 
sal in LElbaz etHTl (|2007l ) extends to densities that are 
a factor of ~ 5 — 6 higher than their lowest densities. 
This would be encompassed by our two lowest density 
bins. Any reversal would therefore be observable in our 
dataset even with our broad density bins, but is not. 

Our results show that the epoch at which a reversal in 
the SFR-density relation takes place, if any, is still under 
debate. Above a mass limit of log M/M© > 10.25, we 
find that the SFR-density relation does not reverse at 
z < 1. Whe n studying ga l axies at z ~ 0.85 in a fixed 
mass range, ICooper et al.l (|2010j l find a larger propor- 
tion of red galaxies at high densities compared to low 
densities. Assuming these red colors reflect older stel- 
lar populations rather than dusty SF, their findings lend 
support to the dec lining SFR-density relation found in 
iPatel et all (|2009bl) and in this work. The debate also ex- 
tends to higher redshifts as studi es of environment have 
begun to expand beyond z > 1. iTran et all 1)2010( 1 find 
evidence for a high frac tion of SFGs associa ted with a 
cluster at z ~ 1.6. Also. lHavashi et al.l ()2010D find [Oil] 
derived SFRs for galaxies in a cluster at z ~ 1.5 that 
are si milar to those of ga laxies in the outskirts. In con- 
trast, |GobaFeF^l] (|2010D find a well formed X-ray emit- 
ting cluster with a large r ed elliptical ga l axy p opulation 
at z ~ 2. Also at z ~ 2, iTanaka et all (|2010D find the 
SFRs of galaxies in a proto-cluster to be lower than that 
of galaxies in the field. Meanwhile, wide-field clustering 
studies at 2 < z < 3 find red galaxies to be strongly 
clustered, suggesting that a color-density relation was in 
place at those early t imes with red galaxies do minating 
more massive halos (jQuadri et al.l 120071 I2008D . These 
works point out however that one key unresolved ques- 
tion is the source of the red colors: dusty SF or old stel- 
lar populations. As a followup, Quadri et al. (2011, in 
preparation) find the fraction of quiescent galaxies, se- 
lected with a UVJ diagram similar to ours, to increase 
with density out to z ~ 2. 

6.2. The Origin of the SFR-Density Relation at 
0.6 < z < 0.9 



One of the key questions we set out to answer was how 
the SFRs of individual galaxies were changing in order to 
produce the declining SSFR-Density relation (Figure [8]). 
Either we are seeing (1) a changing mix of quiescent and 
SFGs with increasing density, (2) a decline in the SFRs 
of SFGs with density, or (3) a combination of the two 
effects. 

From Figure [TTJ we infer that the fraction of SFGs de- 
clines from ~ 68% at the lowest densities to 21% at the 
highest densities, a factor of ~ 3 decrease. If the SFRs 
of SFGs were to remain constant across various environ- 
ments, and the SFRs of quiescent galaxies were negligi- 
ble, then the change in the relative proportion of quies- 
cent and SFGs with environment would contribute only 
a factor of ~ 3 towards the decline in the SSFR-density 
relation. This amount is less than the order of magni- 
tude decline found in Figure |S] for the SED and [Oil] 
based SFRs, but close to the factor of ~ 4 decline found 
for the MIPS based SFRs. The change in the fraction 
of quiescent and SFGs with density therefore cannot be 
solely responsible for driving the declining SED and [Oil] 
SSFR-density relations. Indeed, this is what we find in 
Figure [T2l with the SFRs of SFGs declining by factors of 
~ 5-6. The MIPS based SFRs of SFGs also decline, but 
by a much smaller amount. Thus, both (1) and (2) above 
contribute towards the declining SSFR-density relation 
and potentially reflect a diversity of mechanisms that 
are respo nsible for quench i ng SF on different timescales 
(see, e.g.. lTreu et "aL|[200l iMoran et al.|[2007l ). 

Other environmental studies at somewhat lower red- 
shifts also find the SFRs of SFGs to decline at higher 
densities. For example, using a different color-color se- 
lection from our own (i n addition to othe r selection cri- 
teria) to identify SFGs, IWolf et al.l (pOOflh found the SS- 
FRs of SF Gs to decline at higher densities at z ~ 0.17. 
Meanwhile, Poggi a"n~ti et al.l (|2008l ) argue that SF activ- 
ity in SFGs at 0.4 < z < 0.8 is independent of local 
density based on [Oil] equivalent widths. However, a 
straight line fit to the log(SSFR)-lo g(density) relation in 
Figure 7 of iPoggianti et all (|2008h indicates a negative 
slope of A (log SSFR)/A(log density) « -0.5 at a sig- 
nificance of ~ 4cr. We note that our defini tion for SFGs 
differs from that of IPoggianti et ail (120081) . Also , while 
the luminosity-limited sample of IPoggianti et all (2008) 
is limited to regions in the vicinity of clusters, our survey 
samples a broader range of environments. 

The decline in SFRs with increasing density for SFGs 
seen in Figure [12] could suggest a gradual shut down 
in the cold gas supply that fuels SF at 0.6 < z < 0.9. 
This is further supported by the decreasing value for Ay 
with increasing density. As dust and cold gas are tied 
to the SF process, the declining Ay at higher densities 
suggests less fuel for SF: recent observations of HI de- 
ficient sp_indj5Mii_tlw_V^ dust 
disks dCortese et al.ll2010fc |Pe Looze et alll2010D provide 
strong support for this view. Other recent works at low 
rcdshift also point to the slow remov al of gas in order 
to rep ro duce various observation s (e.g.. fWeinmann et all 
200 9ajbl : Ivan der Wei et alll2010f ). 

Finally, we note that our selection of SFGs in this 
work utilizes the color-color space in which they are pre- 
dicted and observed to occupy. Other works that require 
the detection of emission lines or mid-IR flux find that 
such galaxies have relatively little variation in SFR at a 
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fixed mass (intrinsic 1 — a scatter of ~ 0.3 — 0.5 dex) 
(see, e.g., iBrinchmann et all 120041: iNoeske et ail l2007t 
iRodighiero et al.ll2010f ). Meanwhile, the range in SFRs 
(SED and [Oil]) for our UV J -selected SFGs across our 
three lowest density bins (i.e., environments typically 
probed by the field surveys above) is roughly ~ 0.6 dex. 
Thus, our selection of SFGs and interpretation of their 
properties in this work may differ from those of others. 

6.3. SFRs of SFGs from MIPS vs. Other SFR 
Indicators 

The SSFR-density slope for SFGs in Figure [12] for 
MIPS based SFRs is much shallower compared to the 
SED and [Oil] based SFRs. Several recent results could 
explain these findings. For example, while the overall 
SF activity is clearly diminished in higher density re- 
gions, the remaining S F that takes place may occur in 
a more obscured state. iWolf et al.l (|2009ft find that the 
few SFGs at high densities at z ~ 0.17 that remain are 
predominantly spirals of the reddened v ariety (see also, 
Galla zzi et alJl200ll IHaines et al.l[20lH . IBekki fc ComE 
(120101) propose that such red, spiral SFGs in higher den- 
sity environments result from "outside-in truncation" , 
with the remaining SF taking place in the inner regions 
of the galaxy where higher gas densities and metallic- 
ities lead to heavy obscuration. Another explanation 
could be that increased AGN activity in hig h density 
regions at these redshifts ((Martini et al.|[2009t) could be 
responsible for the elevated mid-IR flux. Our LDP spec- 
tra lack the resolution necessary to identify Type 2 AGN 
from emission line diagnostics and therefore such objects 
likely contaminate our sample. We note however that 
less than < 1% of our sample has an X-ray counterpar t 
from the ChaMP point source catalog (|Kim et al.| [2007), 
and ignorin g these objects does not impact our results. 
In addition. iFu et ail (|2010f) used Spitzer IRS spectra to 
show that objects at z ~ 0.7 with MIPS 24 fim fluxes 
less than ~ 1 mJy, as is the case with our entire sample, 
are dominated by SF rather than AGN activity. Finally, 
the shallow decline of SSFR with density for SFGs in 
Figure [TJ] from the MIPS based SFRs may reflect the 
potentially lo nger timescales probed by that SFR indica - 
tor (see, e.g.. ISalim et al.ll2009t iKelson fc Holdenll2010j ). 
High density regions accrete galaxies from lower densi- 
ties where the typical SFRs are much higher. If a SFR 
indicator traces SF over long timescales, the median SFR 
for a sample of galaxies at high densities will be higher 
given that the newly accreted galaxies from low density 
will still register high SFRs. Lending further support 
to this explanation is the finding that MIPS 24 /jm de- 
tected sources at high dens ities at these redshifts exhibit 
strong Balmer absorption (|Kocevski et al.lfeOlOal ). This 
would indicate a substantial population of A-stars of age 
~ 1 Gyr old, a tim esc ale in accord with the fin dings of 
ISalim et al.l ll2009l) and IKelson fe Holdenl tf2010h . 

7. SUMMARY 

We have carried out a wide-field spectroscopic survey 
of galaxies at 0.6 < z < 0.9 that span a range of en- 
vironments, including the low density field, groups, and 
the z ~ 0.83 cluster RX J0152-13. We selected galaxies 
with zab < 23.3 mag for spectroscopy with IMACS on 
Magellan with a low-dispersion prism (LDP). These LDP 



spectra yielded redshifts with a remarkable precision of 
a s /(l + z)~l%. 

We used a mass-limited sample above 
log M/Mq > 10.25 to study the SFRs of galaxies 
as a function of the local galaxy density. Above the 
mass limit, the distribution of stellar masses in different 
density bins in our sample is very similar (Figure II 1 [) . 
Thus the environmental trends reported in this paper 
are not significantly affected by any residual correlation 
with stellar mass. 

Given the systematics associated with any single SFR 
indicator, we used three different SFR indicators in draw- 
ing conclusions about the environments at 0.6 < z < 0.9 
where SF is quenched. We computed SFRs from: (1) 
rest-frame 12-15 /jm luminosities measured with Spitzer 
MIPS 24 /im imaging, (2) SED fitting, and (3) extinc- 
tion corrected [OII]A3727 A luminosities from the LDP 
spectra. 

To distinguish SFGs from quiescent galaxies we used 
a U — V vs. V — J diagram (referred to here as the 
UVJ diagram). This allowed us to determine whether 
the SFR-density relation at 0.6 < z < 0.9 was driven by 

(1) a change in the proportion of quiescent and SFGs or 

(2) a change in the SFRs of SFGs in different environ- 
ments. The UVJ diagram is especially useful in isolating 
dusty, SFGs on the red-sequence, which would otherwise 
be classified as quiescent based on a single optical color 
(Figures @] and H| • 

Our main results are the following: 

1. With three different SFR indicators, we find that 
the median SSFR-density relation at 0.6 < z < 0.9 
declines across the entire range of densities probed: 
from the low-density field, to the cores of groups 
and a cluster (Figure |8|). The declining SSFR- 
density relation holds even when removing galaxies 
near the cluster (Figure [5]). The declining SSFR- 
density relation is therefore general to all environ- 
ments at these redshifts and not driven by galaxies 
located near the rich cluster in our sample. 

2. For fixed mass bins, the SSFR-density shows a clear 
decline with increasing density (Figure I X0[) . con- 
firming that the local environment is indeed a fac- 
tor in impacting the SFRs of galaxies when con- 
trolling for mass. This result also implies that at 
a fixed mass the SFR-density relation at these red- 
shifts follows a similar decline with increasing den- 
sity. 

3. The fraction of quiescent galaxies increases by a 
factor of - 2.5, from - 32% to - 79%, be- 
tween the lowest and highest densities (Figure ITT)) . 
We note that these quiescent fractions are simi- 
la r to the morphological ea rly type fractions seen 
in lvan der Wei et al.l (|2007t ) when adjusting for the 
mass selection. We also find that the SSFRs of 
SFGs decline with density (Figure Q2]). The de- 
clining SFR-density and SSFR-density relation at 
0.6 < z < 0.9 can therefore be attributed to a com- 
bination of (1) a changing mix of quiescent and 
SFGs, and (2) a decline in the SFRs of SFGs with 
increasing density. 

4. In parallel with the declining SFRs of SFGs, the 



18 



median Ay of SFGs also decline at higher densities 
(Figure [l3f)) . This decline in Ay signifies a lack 
of cold gas and therefore the quenching of SF in 
higher density environments. 

In future work, we plan to use HST ACS and WFC3 
imaging to determine how SF activity is distributed 
within SFGs in different environments. Such an anal- 
ysis will lead to insights on whether SF is being turned 
off across the entire galaxy (e.g., via starvation) or only 
in the outer regions (e.g. gas stripping). The answer has 
implications for the mechanisms that are responsible for 
regulating SF in galaxies residing in different environ- 



ments at half the age of the universe. 
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APPENDIX 

STELLAR MASS LIMIT OF THE SAMPLE 

We primarily use a stellar mass limited sample to conduct our study. Figure [14] illustrates how we determine our 
mass limit. The figure shows the stellar masses of galaxies with zab < 23.3 mag as a function of rcdshift for galaxies 
at 0.6 < z < 0.9 and with K s imaging. Using the UVJ diagram (see Figure 2]), we identify quiescent galaxies and 
denote them as red circles in the figure. SFGs are shown as open blue diamonds. Note how SFGs extend to low masses 
where low mass quiescent counterparts are lacking. In order to study a representative stellar mass limited sample of 
quiescent and SFGs the mass limit will then be determined by the redder quiescent galaxies. The solid line shows the 
stellar mass for an old stellar population formed at Zf = 4 and observed to have a magnitude of zab =23.3 mag at 
all redshifts. Note that a similar galaxy with a fainter magnitude would fall below this line. For a quiescent galaxy 
in our sample at zab = 23.3 mag (i.e., our survey magnitude selection limit) the limiting stellar mass at z = 0.9 is 
logM/Af Q > 10.25 (dashed line). This value represents the limiting stellar mass for galaxies in the redshift interval 
studied in this paper, 0.6 < z < 0.9. 

SPECTROSCOPIC COMPLETENESS 

In making various measurements in this paper, we have taken care to account for variations in the spectroscopic 
completeness of our survey. The completeness here is defined as the ratio of the number of galaxies with a high 
quality LDP redshift to the number of galaxies in the z-band selection catalog. The reciprocal of the completeness 
fraction is used as a weight when computing quantities such as the local galaxy density, median SSFRs, etc. Figure [TBI 
shows the spectroscopic completeness as a function of the observed R — z color and z-band magnitude. Panel (a) 
shows the completeness within an ellipse centered on RAo = 28.1814 deg and Deco = —13.9998 deg, with axis lengths 
ora = 0.174 deg and &Dcc = 0.132 deg. Panel (b) shows the completeness in the outer regions with the same RAo 
and Deco but with <xra = 0.29 deg and &Dcc = 0.22 deg. In general, the inner region has a slightly higher level of 
completeness. Overplotted on each completeness map are the colors and magnitudes of objects within the given region 
and with redshifts in the range 0.6 < z < 0.9. The bulk of our sample lies in regions of completeness space that are 
relatively smooth. 
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Figure 14. Stellar mass vs. redshift for galaxies at 0.6 < z < 0.9 with z^B < 23.3 mag and with K a imaging. J7 VJ-selected quiescent 
galaxies are indicated by the solid red circles, while SFGs are shown as open blue diamonds. The solid black line indicates the stellar 
mass of a galaxy that has undergone an instantaneous burst at Zf = 4 and is observed to have zab = 23.3 mag at all redshifts. This 
model indicates that our sample can be made representative for young SFGs and old quiescent galaxies above a limiting stellar mass of 
logM/Mo > 10.25 (dashed black line). Note that the mass completeness limit for SFGs (blue points) is much lower. 



cn 
O 

M 
I 

(Z 



(a) R<15': 81.3% 




(b) R>15': 69.5% 



I 







1, 






.■'.•■'VS.':" 
























~ • * 



18 



19 



20 
z 



AB 



21 22 23 18 19 20 21 22 23 
(mog) z AB (mog) 



Figure 15. Spectroscopic completeness as a function of observed R — z color and z-band magnitude for (a) the inner R < 15' and 
(b) the outer R > 15' regions of the z-band imaging (see text for precise boundaries of each region). The overall completeness for each 
region is indicated in the top left of each panel. A color bar corresponding to the completeness fraction is shown in the top left. We 
note that the R — z magnitudes represent Vega— AB mag as the raw data were zeropointed using the system that was native to each filter 
(Rab = Kvcga + 0.185 mag). Each panel also shows the R — z color and z-band magnitudes of galaxies in the given region at 0.6 < z < 0.9, 
giving a sense for where most of our sample lies in these completeness maps. 
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